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Abstract 

While metamaterials offer the potential to realize Tamm mediums which represent vacuous spacetime 
subjected to gravitational fields, practical formulations for suitable metamaterials have not hitherto been 
developed. This matter is addressed by establishing a metamaterial formulation for the Tamm medium 
representing Schwarzschild-(anti-)de Sitter spacetime. Our formulation is remarkably simple and does 
not involve a complex nanostructure of the type that is often associated with metamaterials. Instead it 
is based on the homogenization of only isotropic dielectric and isotropic magnetic component mediums, 
which are distributed randomly as oriented spheroidal particles. Using the inverse Bruggeman homoge- 
nization formalism, we demonstrated that a wide range of constitutive parameter values for the Tamm 
medium may be accessed through varying the particle shape, volume fraction or relative permittivity 
and relative permeability of the component mediums. The presented formulation is appropriate for the 
regions of spacetime which lie outside the event horizon for Schwarzschild spacetime and inside the event 
horizon for de Sitter spacetime; there are no such restrictions for anti-de Sitter spacetime. 

PACS numbers: 04.40.Nr, 98.80.-k, 41.20.Jb 

Keywords: metamaterials, Bruggeman homogenization formalism, Tamm medium, Schwarzschild- (anti-)de 
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1 Introduction 

Metamaterials constitute a notable class of engineered materials which offer opportunities for realizing such 
exotic phcnomcnons as negative refraction and cloaking [1]. Furthermore, metamaterials offer unique op- 
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portunities to investigate general relativistic scenarios [2]. This arises from the formal analogy that ex- 
ists between light propagation in vacuous curved spacetime and propagation in a certain nonhomogeneous 
anisotropic or bianisotropic medium, known as a Tamm medium [3, 4, 5]. Lately, theoretical metamatcrial- 
based analogs of black holes [6], de Sitter spacetime [7, 8], strings [9] including cosmic strings [10], and 
wormholes [11], for examples, have been proposed. Crucially, these spacetimes are amenable to representa- 
tion by mctamaterials because their metrics are time- independent. 

While metamaterials may in principle be exploited to construct analogs of curved spacetime, concrete 
details of how this can be achieved in practice are conspicuously absent from the literature. A noteworthy 
exception is a recent description of a metamaterial representation of an artificial electromagnetic black 
hole [12], based upon the homogenization of simple components. However, the two-dimensional black hole 
considered is not astrophysical. In the following sections we develop a metamaterial formulation for a rather 
more complex curved spacetime scenario, which includes the Schwarzschild black hole and (anti-)de Sitter 
spacetime as specializations. Our formulation — which is remarkable for its simplicity — relies on the 
homogenization of isotropic dielectric and isotropic magnetic component mediums, which are distributed 
randomly as oriented spheroidal particles. 

As regards notation, 3- vectors are underlined with the " symbol denoting a unit vector, whereas 3x3 
dyadics are double underlined with / being the identity. The speed of light in vacuum in the absence of a 
gravitational field is Co = l/^e^, where eo = 8.854 x 10"-^^ F m ~^ and /Xo = 47r x 10~^^ H m~^. 



2 Tamm medium for Schwarzschild- (anti-)de Sitter spacetime 

Static Schwarzschild-(anti-)de Sitter spacetime is conventionally represented by the line element [13, 14, 15] 

dS^ = (l-F^ dp - dr^ - f'^idO^ + sin^Od^'^), (1) 

expressed in spherical coordinates with the adopted signature — , — , — ). Herein the function 

~ 2GM Af2 

with M > being the mass of the corresponding Schwarzschild black hole, G the gravitational constant, and 
A the cosmological constant. The spacetime is called Schwarzschild-de Sitter spacetime if the cosmological 
constant A is positive; on the other hand, Schwarzschild-anti-de Sitter spacetime is characterized by A < 0. 
[16]. There are two noteworthy specializations, namely: 

(a) the A = regime corresponds to Schwarzschild spacetime, and 

(b) the M = regime corresponds to de Sitter or anti-de Sitter spacetime, according to whether A is 
positive or negative, respectively. 
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In terms of Cartesian coordinates x = fsin0cos(/), y = f sm0 sin cj), and z = fcosO, the line element (1) is 

1/3 



represented by the metric ga/s as'^ 
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For our purposes here, it is more convenient to work with a diagonal metric. Therefore, we implement the 
(spatial) coordinate transformation represented by 



t x y z)^ = M- (f X y 2)'^. 



(4) 



with the change of basis matrix 
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Since M is an orthogonal matrix, we have f = r = \J x^ + y"^ + z'^. With respect to this new coordinate 
system, the metric has the diagonal form 



fl-F 



V 



\ 

F - 1 
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(6) 



Following the noncovariant approach pioneered by Tamm [3, 4, 5], the electromagnetic response of vacuum 
in curved spacetime represented by the metric (6) may be described by the constitutive relations of an 
equivalent, instantaneously responding, medium (known as a Tamm medium) per [19] 



(7) 



wherein SI units are implemented. Here, 7 is the 3x3 dyadic equivalent of the metric [jab ] with components 



Jab 



goo ' 



(8) 



'Roman indexes take the values 1, 2 and 3; while Greek indexes take the values 0, 1, 2, and 3. 
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That is, the dyadic 7 has the uniaxial form 

2 = diag(l,7,7) , (9) 

where the scalar 7 = 1/ (^1 — F^. 

In view of our quest to construct a Tamm medium with relative permittivity dyadic 7 and relative 
permeability dyadic 7, it is of interest to consider the spatial dependency of the parameter 7. In the case 

of Schwarzschild spacetime, 7 with respect to r exhibits a singularity at r = 2GM/cg, which corresponds 
to an event horizon. In the case of (anti-)de Sitter spacetime, 7 with respect to r exhibits singularities 
at r = ±Co-\/3/A. Since r > 0, this corresponds to a single event horizon for de Sitter spacetime at 
r = Co-y/3/A and no event horizon for anti-dc Sitter spacetime. Furthermore, 7 > outside the event horizon 
for Schwarzschild spacetime but inside the event horizon for dc Sitter spacetime; and 7 > everywhere 
for anti-de Sitter spacetime. These features arc illustrated in Fig. 1, wherein 7 is plotted versus r for 
Schwarzschild and (anti-)dc Sitter spacetimes for M = |A| = 1.3 using the normalizations Cq = G = 1. 

The situation for Schwarzschild- (anti-)dc Sitter spacetime is rather more complicated, as is shown in 
Fig. 2 where 7 is plotted versus r for Schwarzschild-(anti-)de Sitter spacetimes for M = |A| = 0.3 using the 
normalizations Cq = G = 1. Now, for r > 0, 7 with respect to r exhibits two singularities for Scliwarzschild-de 
Sitter spacetime (and a third one for r < 0), while there is also one singularity for Schwarzschild-anti-de 
Sitter spacetime. Also, 7 > between the two singularities for Schwarzschild-de Sitter spacetime and outside 
the singularity for Schwarzschild-anti-de Sitter spacetime. 

The distinction between the 7 > and 7 < regimes is of relevance when one considers the phenomenon 
of electromagnetic planewave propagation with negative phase velocity, which is pertinent to metamaterials 
which support negative refraction [17]. In fact, negative phase velocity can arise in the 7 < regime but 
not in the 7 > regime [18, 19]. The homogenization procedure presented in Sec. 3 is suitable for the 7 > 
regime only. 

Finally in this section, we comment on a practical aspect of realizing the Tamm medium. In practice, the 
nonhomogeneous nature of the Tamm medium could be catered for by subdividing the space of interest into 
local neighbourhoods which are sufBciently small to be considered approximately homogeneous. The inverse 
homogenization procedures developed in Sec. 3 would then be applied locally. This piecewise homogeneous 
approach for the Tamm medium is documented in detail elsewhere [20]. 



3 Inverse Bruggeman formalism 

The Tamm medium described by the constitutive relations (7) is a uniaxial dielectric-magnetic medium with 
identical relative permittivity and relative permeability dyadics. In order to construct such a medium, we 
turn to homogenization using as our basis the well-established Bruggeman formalism [21, 22]. 

Let us consider the homogenization of four component mediums, labelled a, b, c and d. Components 
a and b are isotropic dielectric mediums with relative permittivities and e;,, and relative permeabilities 
fJ-a = fJ-b — ^- Components c and d are isotropic magnetic mediums with relative permeabilities fic and 
fid, and relative permittivities ec = ed = 1- The four component mediums are assumed to be randomly 
distributed, with respective volume fractions /a, ft, fc, fd G (Oj 1) with fd = ^ — fa — fb — fc- Each 
component medium is composed of spheroidal particles which are small compared to the electromagnetic 
wavelengths under consideration. The axis of these spheroids for all four component mediums is taken to be 
aligned with the symmetry axis of 7, namely the x axis. Thus, the surface of each spheroid relative to its 
centre is prescribed by the vector 

r,=pig,-r, (10) 

where the positive-definite shape dyadic 

£^ = diag(l,f/,,C/f), {i^a,b,c,d), (11) 
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the radial unit vector is f, and is a linear measure of size. Spheroids characterized by J/^ > 1 are prolate 
whereas oblate spheroids are characterized by C/^ < 1, and the degenerate spherical case corresponds to 
Ut = l. 

The Brue; geman formalism provides estimates of the relative permittivity dyadic = diag (e^^., esn csr) 
and the relative permeability dyadic /i = diag (/i^^, /is^j /^Sr) of the homogenized composite medium 
(HCM), as follows. Let us introduce the polarizability density dyadics 



iBr 



iBr 



{£ = a, b, c, d) . 



The depolarization dyadics iZ'^'^ herein are given by [23, 24] 

m = diag {or , DT ,DT), {£ = a, b, c,d; m = e, p) , 



where the components 
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and the dimensionless parameter ai = U^mg^/iriBr- Parenthetically, these depolarization dyadics are only 
defined for the at > regime [23] — which corresponds the 7 > regime for the Tamm medium. According 
to the Bruggeman formalism, the constitutive parameters of the HCM are related to those of the component 
mediums by the dyadic equations [21, 22] 

A' =0 ] 

(17) 



where 



fa a: 



(771 = e,^). 



(18) 



In fact, as the diagonal dyadics -A™ have the form diag (A™^, A™, A™), {m = e, /i), the dyadic equations 
(17) contain only four independent scalar equations, which are coupled via the constitutive parameters for 
the HCM. 

Conventionally, homogenization formalisms are applied in the forward sense, wherein the constitutive 
parameters of the HCM are estimated from a knowledge of the constitutive parameters of the component 
mediums. However, since our aim here is to find values of ea,6, /^c,d, fa,b,c and Ua.b,c,d such that the 
corresponding HCM coincides with the Tamm medium specified by the constitutive relations (7), we apply 
the Bruggeman formalism in its inverse sense. While formal expressions of the inverse Bruggeman formalism 
have been developed [25], these formal expressions can be ill-defined [26].^ In practice, it is more convenient 
to exploit direct numerical methods in order to implement the inverse formalism [28]. 

The following three distinct applications of the inverse Bruggeman formalism are considered, and illus- 
trated using numerical examples in the next section. In each application, there are four scalar parameters 
to be determined. 



''We note that certain constitutive parameter regimes have been found to be problematic for the inverse Bruggeman homog- 
enization formalism [27], but these regimes do not overlap with the regimes considered here. 
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(i) Assuming that the relative permittivities ea,b and relative permeabilities /^c,ci are known, and that 
Ua = Ub — Uc — Ud = U, we determine the common shape parameter U and the volume fractions fa, 
ft and fc- 

(ii) Assuming that the relative permittivities ea,b and relative permeabilities Hc,d are known, and that the 
volume fractions fa.b,c are fixed, we determine the shape parameters Ua, Ub, Uc and Ud- 

(iii) Assuming that the shape parameters Ua,b.c.d and the volume fractions fa,b,c arc fixed, wc determine 
the relative permittivities ea,6 and relative permeabilities fic.d- 

As a representative example, let us concentrate on the numerical implementation of application (i) — 
the numerical implementations for applications (ii) and (iii) are analogous. A modified Newton-Raphson 
technique [29, 30] may be implemented to extract the volume fractions fa,b.c and common shape factor U 

from cqs. (17). Using this technique, the solutions at step fc + 1, namely ^U''''^^\ fa''^^\ flj^'^^\ fc^^^^^, 
are derived from those at step fc, namely : fa \ fb \ fc I , via the recursive scheme 
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(19) 



wherein the components of the dyadics A'^'^ are expressed as functions of the unknown parameters {[/, fa, fb, fc}- 

For convergence of the scheme (19), it is vital that the initial estimates arc chosen 

to be sufficiently close to the true solution. The forward Bruggcman formalism can be utilized in order to 
generate suitable initial estimates, as we now outline. Let {e^^, tsr, fiBr' I^Br} denote the forward Bruggeman 
estimates of the HCM's relative permittivity and relative permeability parameters, computed for physically 
reasonable ranges of the parameters U and fa,b,c, namely U € (C/+, U~) and fa,b,c G (^fab c fab cj ■ Then: 

(1) Fix fa - (/- + /+) /2, fb = if,- + /+) /2, and = (/- + /+) /2. For all values of [/ e {U-,U+), 
find the value W for which the quantity 



1)^ 



^Br 



(fiBr - 1) + 



1 



(20) 



is minimized. 



(2) Fix U = U\ fb = {fb + fb) /2, and f, = (/- + /+) /2. For all values of fa e {fa, ft), find the value 
f\ for which A is minimized. 

(3) Fix U = f/t, fa = fl and f, = (/- + /+) /2. For aU values of fb £ {fbJb)^ ^iid the value fl for 
which A is minimized. 

(4) Fix U = U\ fa = fl, and fb = fl For ah values of fc £ ifcJc)^ find the value /] for which A is 
minimized. 
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The steps (l)-(4) are repeated, using /t, and as the fixed values of fa,b,c in step (i), and as the 
fixed values of /f,_c in step (ii), and as the fixed value of fc in step (iii), until A becomes sufficiently small. 

In our numerical studies, we found that when A < 0.01, the values of U\ f^, fl and fl provide suitable 
initial estimates for the modified Newton- Raphson scheme (19). In fact, we found that this technique 
of iteratively scanning the space of possible solutions could itself be used to find the inverse Bruggeman 
solutions, in some cases with a faster rate of convergence than the Newton-Raphson method and in a 
manner that is less sensitive to the initial estimates. 

4 Numerical illustrations 

We present numerical illustrations of the inverse homogenization applications (i)-(iii) described in Sec. 3. 
For each illustration, the parameter estimates for the component mediums are calculated as functions of 7. 
The ranges 7 > 1 and < 7 < 1 are considered — the 7 > 1 range corresponds to Schwarzschild-dc Sitter 
spacetimc whereas the < 7 < 1 corresponds to both Schwarzschild-de Sitter and Schwarzschild-anti-de 
Sitter spacetime, depending upon the magnitude of A. For all the numerical results presented in Figs. 3-5, 
the degree of convergence of the numerical schemes was < 1%, and for most of the plotted points this value 
was < 0.1%. 

Let us begin with application (i). The common shape parameter U and the volume fractions fa.b,c are 
plotted versus 7 in Fig. 3. For 0.8 < 7 < 0.9 we chose Ca = 2, £{, = 0.1, = 1.8 and /id = 0.2, while for 
1 < 7 < 2.9 we chose Ea = 8, = 0.3, /ic = 7.4 and /i^ = 0.4. The common shape parameter increases 
sharply as 7 increases, for both 0.8 < 7 < 0.9 and 1 < 7 < 2.9 ranges, while no particularly noteworthy 
trend is obvious from the volume fraction plots. 

In Fig. 4, the shape parameters Ua.b,c,d are presented as functions of 7 for the application (ii). For 
0.795 < 7 < 0.805 we chose = 2, = 0.1, /ic = 1-8 and fid = 0.2 (as we chose for Fig. 3 for 0.8 < 7 < 0.9) 
and fa.b = 0.21, fc = 0.3 , while for 1.9 < 7 < 2.2 we chose Sa ~ 8, ei, = 0.3, fic = 7.4 and /i^ = 0.4 (as 
we chose for Fig. 3 for 1 < 7 < 2.9) and fa.b,c ~ 0.25. We see that Ua,c decrease sharply as 7 increases for 
0.795 < 7 < 0.805, but increase sharply as 7 increases for 1.9 < 7 < 2.2; the opposite trend is exhibited by 

f/M- 

Finally, for application (iii), plots of the relative permittivities ea,6 and relative permeabilities fic,d against 
7 are provided in Fig. 5. For 0.7 < 7 < 1 we chose the common shape parameter U = Ua,b,c,d = 0.01 and 
the common volume fraction fa,b.c = 0.25, while for 1.1 < 7 < 2.9 we chose the common shape parameter 
U = Ua^b.c.d ~ 5 and the common volume fraction /a,h,c — 0.25. For both 0.7 < 7 < 1 and 1.1 < 7 < 2.9 we 
found that the numerical schemes provide values of and /i^ which arc almost the same, and values of ei, 
and fid which are almost the same. 

We close this section by remarking that the relative permittivities and relative permeabilities featured 
in Figs. 3-5 are neither unreasonably high nor unreasonably low from a physical perspective. Indeed, 
as regards the smallest values of ta,b,c,d and fia,b,c,d featured in Figs. 3-5, we note that materials with 
relative permittivities and relative permeabilities close to zero are currently subjects of intense investigation 
[31, 32, 33]. 

5 Closing remarks 

By the homogenization of simple arrangements of four simple component mediums, the Tamm medium for 
Schwarzschild-(anti-)de Sitter spacetime may be constructed. The simplicity of this construction is espe- 
cially noteworthy in view of the complex nanostructures which characterize many metamaterials designed 
for cloaking or negative refraction applications [1], for example. A wide range of constitutive parameter 
values for the Tamm medium can be accessed through varying the particle shape, volume fraction or relative 
permittivity and relative permeability of the component mediums. Thus, the inverse homogenization formu- 
lation delivers a practical strategy for designing an experimental analog for Schwarzschild-(anti-)de Sitter 
spacetime. 
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While the inverse homogenization scenario chosen for presentation in Sec. 3 was based on four component 
mediums, this formulation is not unique. Indeed, fewer components mediums may be used. For example, the 
Tamm medium for Schwarzschild-(anti-)dc Sitter spacctimc could similarly be conceptualized as arising from 
the homogenization of two component mediums, both of which being isotropic dielectric-magnetic mediums 
distributed as oriented spheroids. Alternatively, two component mediums could used which were both uniax- 
ial dielectric-magnetic mediums (with parallel symmetry axes) distributed as spherical particles. However, 
the four-component formulation presented herein incorporates very simple components, and accordingly 
offers a very large degree of freedom in choosing constitutive parameters. 

Lastly, let us comment on the restriction of the inverse homogenization approach to those regions of 
spacetime which lie outside the event horizon for Schwarzschild spacetime and inside the event horizon for 
de Sitter spacetime (there being no such restrictions for anti-de Sitter spacetime). Equivalently, only Tamm 
mediums with positive-definite relative permittivity and relative permeability dyadics can be realized using 
this approach. The constitutive dyadics for the Tamm medium representing Schwarzschild-(anti-)de Sitter 
spacetime cannot possibly be negative definite but, as illustrated in Figs. 1 and 2, the indefinite scenario can 
arise (except in the case of anti-de Sitter spacetime). In principle, metamaterials with indefinite constitu- 
tive dyadics could be harnessed to construct a Tamm medium for the 7 < regime. Indeed, experimental 
studies of anisotropic dielectric [34] and anisotropic magnetic [35, 36] metamaterials with indefinite consti- 
tutive dyadics have already been reported. However, there are likely to be major practical difficulties in 
using such metamaterials to construct a 7 < Tamm medium: First, not only are the relative permittivity 
dyadic and the relative permeability dyadic both required to be indefinite (since the Tamm medium is an 
anisotropic dielectric- magnetic medium), but these two dyadics are required to be identical. Second, the 
metamaterials should be approximately nondissipative in order to faithfully represent the Tamm medium, 
but dissipation has proved to be a severe hindrance for many metamaterials. Finally, we add that con- 
ventional depolarization-dyadic-based homogenization techniques (as typified by the Bruggeman formalism) 
cannot be used to realize such metamaterials as nondissipative HCMs, since the indefinite nature of the 
constitutive dyadics renders the depolarization dyadics undefined [23]. This difficulty could be sidestepped 
by the incorporation of a small amount of dissipation [37], but at the cost of compromising the extent to 
which the resulting HCM represents the desired Tamm medium. 
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support of his research activities. 

References 

[1] T. J. Cui, D. R. Smith, and R. Liu, Eds., Metamaterials: Theory, Design, and Applications (Springer, 
New york, 2010). 

[2] I. I. Smolyaninov, J. Opt. 13, 024004 (2011). 

[3] G. V. Skrotskn, Soviet Phys.-Dokl. 2, 226 (1957). 

[4] J. Plebanski, Phys. Rev. 118, 1396 (1960). 

[5] W. Schleich and M. O. Scully, in New Trends in Atomic Physics, Eds. G. Grynbcrg and R. Stora 
(Elsevier, Amsterdam, 1984), pp.995-1124. 

[6] I. I. Smolyaninov , New J. Phys. 5, 147 (2003). 

[7] M. Li, R.-X. Miao, and Y. Pang, Phys. Lett. B 689, 55 (2010). 



8 



[8] M. Li, R.-X. Miao, and Y. Pang, Opt. Express 18, 9026 (2010). 

[9] R.-X. Miao, R. Zheng, and M. Li, Phys. Lett. B 696, 550 (2011). 
[10] T. G. Mackay and A. Lakhtakia, Phys. Lett. A 374, 2305 (2010). 
[11] A. Greenleaf, Y. Kurylev, M. Lassas, and G. Uhlniann, SIAM Rev. 51, 3 (2009). 
[12] W. Lu, J. Jin, Z. Lin, and H. Chen, J. Appl. Phys. 108, 064517 (2010). 
[13] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2738 (1977). 
[14] J. Podolsky, Gen. Rcl. Grav. 31, 1703 (1999). 
[15] V. Cardoso and J. P. S. Lemos, Phys. Rev. D 64, 084017 (2001). 

[16] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Herlt, Exact Solutions of Einstein's 
Field Equations, 2nd ed. (Cambridge University Press, Cambridge, UK, 2003), chapt. 8. 

[17] T. G. Mackay and A. Lakhtakia, Phys. Rev. B 79, 235121 (2009). 

[18] T. G. Mackay, S. Setiawan, and A. Lakhtakia, Eur. Phys. J. C 41 (Supplement 1), 1 (2005). 
[19] T. G. Mackay, A. Lakhtakia, and S. Setiawan, Europhys. Lett. 71, 925 (2005). 
[20] T. G. Mackay, A. Lakhtakia, and S. Setiawan, New J. Phys. 7, 75 (2005). 

[21] W. S. Weiglhofer, A. Lakhtakia, and B. Michel, Microw. Opt. Technol. Lett. 15, 263 (1997). Corrections: 
22, 221 (1999). 

[22] T. G. Mackay and A. Lakhtakia, Electromagnetic Anisotropy and Bianisotropy: A Field Guide (Word 
Scientific, Singapore, 2010), chap. 6. 

[23] B. Michel, Int. J. Appl. Electromagn. Mech. 8, 219 (1997). 

[24] B. Michel and W. S. Weiglhofer, Arch. Elektron Ubcrtrag. 51, 219 (1997). Corrections: 52, 310 (1998). 
[25] W. S. Weiglhofer, Microw. Opt. Technol. Lett. 28, 421 (2001). 
[26] E. Cherkaev, Inverse Problems 17, 1203 (2001). 

[27] S. S. Jamaian and T. G. Mackay, J. Nanophotonics 4, 043510 (2010). 

[28] T. G. Mackay and A. Lakhtakia, J. Nanophotonics 4, 041535 (2010). 

[29] P. A. Stark, Introduction to Numerical Methods (Macmillan, New York, 1970). 

[30] R. D. Kampia and A. Lakhtakia, J. Phys. D: Appl. Phys. 25, 1390 (1992). 

[31] A. Alu, M. Silveirinha, A. Salandrino, and N. Engheta, Phys. Rev. B 75, 155410 (2007). 

[32] G. Lovat, P. Burghignoli, F. Capolino, and D. R. Jackson, lET Microw. Antennas Propagat. 1. 177 
(2007). 

[33] M. N. Cia, M. Beruete, I. Campillo, and M. Sorolla, Phys Rev. B (accepted for publication). 
[34] H. Lee, Z. Liu, Y. Xiong, C. Sun, and X. Zhang, Opt. Express 15, 15886 (2007). 

[35] A. B. Kozyrev, C. Qin, I. V. Shadrivov, Y. S. Kivshar, I. L. Chuang, and D. W. van dcr Wcidc. Opt. 
Express 15, 11714 (2007). 

[36] J. Sun, L. Kang, R. Wang, L. Liu, L. Sun and J. Zhou, New J. Phys. 12, 083020 (2010). 

[37] T. G. Mackay, A. Lakhtakia, and R. A. Depine, Microw. Opt. Technol. Lett. 48, 363 (2006). 



9 




1 2 3 4 5 

r 

Figure 1: The parameter 7 plotted versus r for Schwarzschild (solid, red), de Sitter (dashed, green) and 
anti-de Sitter (broken dashed, blue) spacetimes. Here M <E {0,1.3}, A G {0,±1.3} with the normalization 
G = Co = 1. 
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Figure 3: The common shape parameter U (thick sohd, red) and volume fractions fa (dashed, green), ft 
(broken dashed, blue) and fc (thin solid, blue) plotted versus 7. The relative permittivities = 2, £{, = 0.1 
and relative permeabilities yUc = 1-8, = 0.2 for 0.8 < 7 < 0.9 (upper); and = 8, = 0.3, ~ 7.4, 
fid = 0.4 for 1 < 7 < 2.9 (lower). 
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Figure 4: The shape parameters Ua (thick soHd, red), Ub (dashed, green), Uc (broken dashed, blue) and Ud 
(thin sohd, blue) plotted versus 7. The relative permittivities = 2, Cb = 0.1 and relative permeabilities 
yUc = 1-8, /id = 0.2 and volume fractions = 0.21, /j, = 0.21 and fc = 0.3 for 0.795 < 7 < 0.805 (upper); 
and ea 8, eb = 0.3, = 7.4, /i^ ^ 0.4 and fa,b,c4 = 0.25 for 1.9 < 7 < 2.2 (lower). 
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Figure 5: The relative permittivities (thiek solid, red), (dashed, green) and relative permeabilities fic 
(broken dashed, blue), fid (thin solid, blue) plotted versus 7. The shape parameter U — Ua.b,c,d — 0.01 
and volume fractions fa,b.c = 0.25 for 0.7 < 7 < 1 (upper); and U ~ Ua,b.c.d — 5 and fa,b.c = 0.25 for 
1.1 < 7 < 2.9 (lower). 
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